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Adiabatic approximation and fluctuations in exciton-polariton condensates 
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We study the relation between the models commonly used to describe the dynamics of nonreso- 
nantly pumped exciton-polariton condensates, namely the ones described by the complex Ginzburg- 
Landau equation, and by the open-dissipative Gross-Pitaevskii equation including a separate equa¬ 
tion for the reservoir density. In particular, we focus on the validity of the adiabatic approximation 
that allows to reduce the coupled condensate-reservoir dynamics to a single partial differential equa¬ 
tion. We find that the adiabatic approximation consists of three independent analytical conditions 
that have to be fulfilled simultaneously. By investigating stochastic versions of the two correspond¬ 
ing models, we verify that the breakdown of the adiabatic approximation can lead to discrepancies 
in correlation lengths and distributions of fluctuations. Additionally, we consider the phase diffusion 
and number fluctuations of a condensate in a box, and show that self-consistent description requires 
treatment beyond the typical Bogoliubov approximation. 
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I. INTRODUCTION 

Exciton-polaritons are coherent superpositions of 
quantum well excitons and a microcavity photons, re¬ 
sulting from strong coupling of the two modes at reso¬ 
nance [1-3]. The mixed nature of these quasiparticles is 
attractive from the point of view of fundamental and ap¬ 
plied research in many ways. The matter component pro¬ 
vides strong interactions, while the photonic component 
yields very light effective mass and allows for straightfor¬ 
ward detection. Condensation and supefluidity of polari- 
tons or polariton lasing has been demonstrated in many 
laboratories, even at room temperature [4-11]. 

In the case of nonresonant pumping, polaritons can be 
created either using a beam at frequency above the po¬ 
lariton resonance or using electrical carrier injection [12]. 
Free carriers and high-energy excitons undergo energy re¬ 
laxation towards the polariton ground state, where they 
can condense [7]. Theoretical description of this com¬ 
plicated process, involving scattering with phonons and 
other polaritons is a formidable task, and several ap¬ 
proaches have been proposed in the past with various 
approximations involved [13-18]. Among these, phe¬ 
nomenological models based on various generalizations 
of the Gross-Pitaevskii equation [19-21] have been par¬ 
ticularly useful thanks to their simplicity and limited 
number of external parameters. The most commonly 
used are the complex Ginzburg-Landau equation (CGLE, 
also termed the generalized GP equation in some works), 
with a single equation for the condensate evolution, and 
the open-dissipative Gross-Pitaevskii equation (ODGPE, 
other names are also used in the literature) [20] , which at¬ 
tempts to describe the full system dynamics using a pair 
of coupled condensate-reservoir equations. It has been 
pointed out in several places [22-24] that the ODGPE 
model can be reduced to the GGLE model under the 
adiabatic (quickly responding reservoir) approximation, 
but the validity of this approximation has not been in¬ 
vestigated in detail. 

In this paper, we investigate systematically the rela¬ 


tion between the the GGLE the ODGPE models. We 
establish precisely conditions under which the reduction 
to the GGLE model is justihed; contrary to the common 
belief, we show that the fast reservoir relaxation alone 
is not a sufficient condition. We show that adiabatic- 
ity requires three independent analytical conditions to 
be fulfilled simultaneously. Additionally, the condensate 
must remain close to the steady state, since large fluctu¬ 
ations may lead to complete breakdown of the correspon¬ 
dence between the models. Such large fluctuations occur 
in particular close to the condensation/stability limits of 
the condensate phase diagram. 

In the second part of the paper, we investigate how the 
breakdown of the adiabaticity influences the steady state 
solutions of the corresponding stochastic GGLE/ODGPE 
models. Recently, fluctuations of nonequilibrium quan¬ 
tum fluids became a very active area of research. Spa¬ 
tial [25, 26] and temporal [27, 28] correlations have been 
investigated in the small fluctuations regime, and the 
critical scaling properties have been established using the 
renormalization group [28-30]. Dynamics of the polari¬ 
ton condensation phase transition have been shown to 
display similarities with the Kibble-Zurek theory of uni¬ 
versal dynamics [31, 32]. Here, we show how the spec¬ 
trum of fluctuations and spatial correlations are modified 
in the non-adiabatic regime as well as in the large fluctu¬ 
ations regime. In particular, we demonstrate the appear¬ 
ance of dark soliton-like structures and chaotic bistable 
steady states close to the limits of condensate stability in 
parameter space in the ODGPE model. 

In addition, we consider a model of a condensate 
in a box, in which case condensation may occur de¬ 
spite the absence of true long-range order at low di¬ 
mensions [25, 33, 34]. We show that a self-consistent 
description requires treatment beyond the typical Bo¬ 
goliubov approximation, and demonstrate how the zero- 
momentum singularity of the momentum distribution 
spectrum [25] can be avoided by an appropriate gener¬ 
alization of the Bogoliubov ansatz. This allows for the 
determination of the number fluctuations and condensate 
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phase diffusion equation. 

The paper is organized as follows. In Sec. II we in¬ 
troduce the CGLE and ODGPE models as well as their 
stochastic versions, in particular the stochastic Gross- 
Pitaevskii (SGPE) model being a generalization of the 
GGLE. In Sec. Ill we derive carefully the conditions un¬ 
der which the adiabatic reduction of the ODGPE to the 
GGLE model is justified. In Sec. IV we recall briefly the 
main analytical results concerning the fluctuations and 
spatial correlations of the condensate. We also present 
the analysis of the fluctuations and phase diffusion of 
a condensate confined in a box. In Sec. V we present 
numerical investigation of the properties of stochastic 
steady states in various regimes. Sec. VI concludes the 
paper. 


II. MODELS 

A. Complex Ginzburg-Landau equation 


loq is the oscillation frequency, Pq is the pumping rate, 
7 o is the polariton loss rate, rig is the saturation density, 
and dIEsGPE(x, f) is the complex stochastic noise. 

As in the second part of the manuscript we will fo¬ 
cus on the spatial correlations, we will compare our re¬ 
sults to the above model (2), which was used in [20, 25], 
rather than the standard GGLE (I). These two mod¬ 
els are, however, completely equivalent as long as one 
is interested in a state with small density fluctuations 
around the homogeneous solution. To show this con¬ 
sider the solution for which the density can be written 
as |0(x, = |(/)op -I- Sn, where i5n/|(/)oP -C 1. The 
fraction in Eq. (2) can be expanded in Taylor series, 
and by neglecting second- and higher-order terms in Sn 
we obtain Eq. (1) with the correspondence A = cjq, 
B = hjlml, C = go, D = ^0 - ijl/Po), E = j^/nsPo, 
dWcGLE = dlVsGPE, which yields the steady state den¬ 
sity = ns{Po - 7o)/7o for Pq > 7o- 

For completeness, we note that in [20, 25] the stochas¬ 
tic noise was assumed to be Gaussian with the correla¬ 
tions 


The simplest dissipative model to describe one¬ 
dimensional exciton-polariton condensates is the 
complex-Ginzburg-Landau equation (CGLE) 

id(j) = [A - -b C\ct>f +i{D- E\cl)f)] (/>df-b 
-bdWcGLE, (1) 


where for future reference we additionally included the 
stochastic term dIVcGLE(x, t) which vanishes in the clas¬ 
sical limit. The real parameters of the equation de¬ 
scribe the energy offset (A), the dispersion coefficient 
(B), nonlinear interactions (C), external pumping {D), 
and nonlinear losses {E). In the case of interest the 
parameters B, C, D and E are positive, which guar¬ 
antees the existence of a st able homogeneous solution 
^(x,t) = D/Ee~''^*' with g = A + C\(j)o\'^ 

[35]. Compared to the most general form of the CGLE, 
we assume no diffusive term (as B is real), although 
models including the diffusive term have also been em¬ 
ployed [19, 29, 36]. It is, however, not crucial to our 
considerations. We note that the coefficient A can be 
removed simply by moving to a rotating frame where 
^'(x, t) = (^(x, t)e*^*, provided that the noise is not time- 
correlated. 

A generalized form of the CGLE, or the stochastic 
Gross-Pitaevskii equation (SGPE), has been used in sev¬ 
eral works which investigated spatial correlations of the 
condensate [20, 25] 


*d(/) = 


^0 - -b5o|</>P + 
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2 m, 

Pn 


( 2 ) 


1 + ^ 


■ 7o 


<j)dt + dWsGPE, 


where ^(x, t) is the wave function, mg is the effective 
mass of lower polaritons, go is the interaction coefficient. 


(dtVsGPE(x)dlVgGpE(x')) = 2D00(5(x - x')dt, 
(dWsGPE(x)dlVsGPE(x')) = 0 (3) 

and the amplitude k, 79 , which is an estimate of 
the quantum noise due to the dissipative coupling. Since 
the above white noise dWsGPE has a diverging norm, 
in practice an appropriate UV cutoff must be employed 
eg. through frequency dependence of the amplification 
term [25]. Alternatively, one may model the process on 
a discretized mesh with lattice constant Ax, where the 
Dirac delta is replaced as ^(x — x') —>• S^.x'/. We 
note that an alternative form of stochatic noise was de¬ 
rived from a quantum model of condensation in a cavity 
interacting with two-level emitters [37]. 


B. Open-dissipative Gross-Pitaevskii equation 


A more realistic model of the polariton condensate, 
called the open-dissipative Gross-Pitaevskii equation 
(ODGPE), includes coupling of the condensate wave 
function to the reservoir described by a density field 
nR(x,t) [14,20] 


idiS = 


2m* 


■ffcjV'P +9RnR+ 


+ 2 - "fc) 


ijjdt -b dW, 


dnn 

dt 


=p - {irAR\'^?) Rr, 


(4) 


(5) 


where the complex stochastic noise can be obtained 
within the truncated Wigner approximation [14] 


(dW(x)dW*(x')) = ^^(RRR + lc)Sx.x'. (6) 
(dlT(x)dW(x')) = 0. (7) 
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where we assumed a constant scattering rate R which 
gives noise correlations analogous as in Eq. (3). Here 
P(x) is the exciton creation rate determined by the 
pumping profile, 7 c and 7 r are the polariton and reser¬ 
voir loss rates, and R is the rate of stimulated scatter¬ 
ing from the reservoir to the condensate, and gc, 5 r are 
the rates of repulsive polariton-polariton and reservoir- 
polariton interactions, respectively. 

In the absence of noise, a spatially uniform solution 
is given by nR(x, t) = n^. Above 

the threshold pumping P > Pth = IcIr/R a stable 
condensate exists with the condensate density |'(/'oP = 
iP/'lc)-hR/R),n% = jclR, and go = ScIV’oP 
We note that formally similar models were also used 
to describe nonlinear effects in semiconductor microcav¬ 
ities at weak coupling [38] and atom lasers [39]. The 
above phenomenological model has been successful in de¬ 
scribing a number of different experimental situations in 
exciton-polariton condensates, although various values of 
parameters have been used in the literature [36, 40-42]. 


III. ADIABATIC APPROXIMATION AND THE 
CORRESPONDENCE BETWEEN THE MODELS 


due to large fluctuations. The full correspondence be¬ 
tween the three models is given by 


. / 27 c 7 c7r\ 


B = 


Wo = Qr 
h 


2ml 


2 m*’ 


f. 7c 

P = go = gc - 9 r-^^ 

Po 2 2PR ’ 


E = 






risPo 2P’ 
-D00 = 70 = 7c/2. 


( 10 ) 


Note that there are only six equations while seven pa¬ 
rameters are present in the open-dissipative model. Con¬ 
sequently, one free parameter can be chosen when cal¬ 
culating seven parameters of Eqs. (4)-(5) corresponding 
to six parameters of Eq. (1) or Eq. (2). In the classical 
limit (no stochastic noise) the last equation is absent and 
there is an additional free parameter in both Eqs. (4)-(5) 
and Eq. (2). 


A. Adiabatic limit 


B. Validity of the adiabatic approximation 


The open-dissipative Gross-Pitaevskii model Eqs. (4)- 
(5) can be reduced to the simpler and more tractable 
CGLE model Eq. (1) under the assumption that the 
reservoir density ^^(x, t) adiabatically follows the change 
of lV^(x, t)j^. With the use of Eq. (5) we can express riR, 
as 


«R(x,f) 


P 

7R-f PlV'(x,t)12' 


( 8 ) 


Note that the above relation is valid also in the case of a 
stationary state nR(x,t) = nR(x), = ^(x)e“*^‘, 

regardless of the adiabatic assumption. However, in this 
work we are interested in dynamical processes and the 
limits of validity of the adiabatic approximation. The 
Eq. (4) now takes the form of a generalized CGLE 


idtp = 
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7R-f PlV’(x,t)12 
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RP 


2 \"fR + R\tlj{x,t)\^ 


■ 7c 


4!^ + dW. (9) 


Again, this equation is equivalent to Eqs. (I) and (2), 
provided that the condensate density j'0(x, t)j^ is close 
to the homogeneous solution \tpo\^ = {P/lc) — (ir/R) 
which allows for expansion of the right hand side in Tay¬ 
lor series to the first order. We emphasize that the expan¬ 
sion is done around the steady state density, and not zero 
density of the condensate. As will be shown in subsec¬ 
tion V C, at low density (pumping slightly above thresh¬ 
old) the results of the two models disagree qualitatively 


The limits of the validity of the adiabatic approxima¬ 
tion can be estimated by comparing the characteristic 
timescales existing in the dynamical system. To this end 
consider Eqs. (4)-(5) and assume that the adiabatic ap¬ 
proximation is fulfilled; the reservoir density tir is able 
to quickly adjust to the condensate density Ji/;]^ which 
changes on a much longer timescale. Equation (5) has a 
simple solution if is treated as a constant 


«R(x,t) 


P 

7r -f Pl'0(x,t)12 




Therefore the timescale of the reservoir is r/j = (y/j -|- 
Pj'0(x, t)j^)“^. For the adiabatic assumption to be con¬ 
sistent, this timescale must be much smaller than the 
timescales of the condensate, which are given by all the 
four terms in Eq. (4). We obtain four conditions that, 
if fulfilled simultaneously, give a sufficient condition for 
the validity of the adiabatic approximation; namely, all 
the terms ^iPriR-yc], ffcjV'Pi and gnUR must be 
much smaller than = jr + R\4’\‘^. The first condi¬ 
tion means that only the low momentum modes of the 
condensate 


k'^ 2m*/ (Htr) (12) 

may be occupied considerably. The other conditions pro¬ 
vide relations between the system parameters, as we show 
below. 

We now make the same assumption that was used to 
compare the different models of the condensate, express¬ 
ing the condensate density as a steady state value plus 
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small fluctuations, |(^(x,t)p = \(j>Q\‘^ + Sn. Under this as¬ 
sumption, the second of the above conditions is always 
fulhlled. The third and fourth conditions take the simple 
forms 


Pth 

0 

1 

(13) 

- 5 -» 

p 

9c 

P ^ 9r 1C 

(14) 

Pth 

R IR 


Equation (13) is always true ii g < R, and otherwise it 
gives an upper limit for the pumping P. Equation (14) 
gives a lower limit for P. Note that the condition (13) is 
independent of the reservoir relaxation rate jr. 

Finally, we note that while the above conditions are 
sufficient, but not necessary, it is unlikely that if one of 
the four terms is large, it could be fully compensated 
by another term in Eq. (4), since the terms depend very 
differently on the fields ur, and ip. We also emphasize 
once again that Eqs. (13) and (14) as well as the exact 
correspondence between the models are only valid if the 
condensate density is close to the steady state value. 


C. Limits of modulational stability 

It is well known that the CGLE (1) displays linear 
modulational instability (also called Benjamin-Feir in¬ 
stability) for BC < 0, which is a special case of the 
Benjamin-Feir-Newell criterion [35, 43]. While the dis¬ 
persion coefficient B is always positive in our model, it 
follows from Eqs. (10) that not all choices of the inter¬ 
action coefficients gc > 0 and gn > 0 correspond to a 
positive value of C. A simple application of Eqs. (10) 
yields the criterion for linear stability 

• ( 15 ) 

^th gc iR 

In other words, the violation of the stability condi¬ 
tion (15) corresponds to effectively attractive interactions 
in the CGLE or SGPE model, C,go < 0. The above in¬ 
stability of the model (4)-(5) was first reported in [20], 
where it was named the “hole-burning effect”, and the 
analytical condition for stability (15) was first derived 
in [44] by investigation of the Bogoliubov-de Gennes ex¬ 
citation spectrum. Alternatively, the same condition can 
be derived from the analysis of the total blueshift origi¬ 
nating from both gc and gn [32]. 


IV. FLUCTUATIONS AND SPATIAL 
CORRELATIONS 


We also derive self-consistently analytical formulas for 
the number fluctuations and phase diffusion of a conden¬ 
sate in a hnite box. In the last section we present nu¬ 
merical results in both the adiabatic and non-adiabatic 
regimes. 


A. Previous analytical results 


In [25] Bogoliubov-de Gennes equations for the small 
fluctuations around the condensate in the form ^(x, t) = 
[00 + ^0(x, f)] were used to determine the momen¬ 
tum distribution of the particles in the SGPE model (2). 
The number of particles in each mode in the steady state 
was calculated as 


n 


SS 

k 




~r~ 


- El 



(16) 


where P = 70 (To ~ 7o)/Tb is the effective damping rate, 
fjL = go 1 00 P > 0 is the energy of interactions between 
particles, and El = ek(ek -I- 2g) with Ck = /2m* is 

the Bogoliubov dispersion. The first term on the right 
hand side corresponds to the condensate, and the second 
term to the non-condensed cloud. 

The formula (16) must be used with caution. It dis¬ 
plays both the IR divergence in the term and the 
UV divergence in the constant term. While the UV diver¬ 
gence can be healed in the case of frequency-dependent 
pumping [25] or discretization of the system on a lattice, 
the IR divergence poses a problem for the interpretation 
of the above formula at the k = 0 point. We show be¬ 
low how this problem can be solved by an appropriate 
generalization of the Bogoliubov approximation. 

Notice also that in the infinite system (“thermody¬ 
namic” limit with U —> 00 at constant pumping) the frac¬ 
tion of particles in the condensate vanishes in the relevant 
one- and two-dimensional cases [25] , invalidating the Bo¬ 
goliubov approximation. However, the long-range corre¬ 
lations are mostly influenced by the phase fluctuations, 
which are much less energetically costly than the am¬ 
plitude fluctuations due to the existence of a Goldstone 
mode of phase twists [33, 34]. The correlations are accu¬ 
rately described using the density-phase generalized Bo¬ 
goliubov approximation 0(x, f) = 

which yields the long-range correlations in one dimen¬ 
sion falling exponentially as [25, 45] 


g^^\x) ^ exp(—x/li), where h 


2noligT 

+ r^) 

(17) 


B. Condensate in a box 


In the rest of the paper, we investigate how the break¬ 
down of the adiabatic approximation influences the fluc¬ 
tuations around the steady state and spatial correlation 
functions. We first recall the relevant results of [25] 
and [45] corresponding to the SGPE (or adiabatic) case. 


In a finite system where k is discretized, the fraction of 
particles in the condensate k = 0 can be large, and we will 
consider this case in this subsection. In the following we 
consider a condensate in a box of length L and periodic 
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boundary conditions with AA: = 27r/L, for which we have 
the correspondence —>• (L/(27r))‘^(5fc_o- The first 

term on the right hand side of (16) can be written simply 
as {2Try\(j)o\‘^d^‘^\k) = |<^opT'^4,o = A^o4,o- To cope 
with the IR divergence due to the singularity at fc = 0 we 
introduce a generalized density-phase Bogoliubov Ansatz 
which takes into account the diffusion of the condensate 
phase (p{t) 

t) = [(j)o + S(j)k=o{t)] t), (18) 


where ^pit) and 5(pk=o(t) are real functions describing the 
phase and amplitude fluctuations of the condensate, and 
5(pkpoi'^,t) is a complex function describing the fluctua¬ 
tions in the out-of-condensate modes. In contrast to the 
standard Bogoliubov approach, we treat the fluctuations 
of the k = 0 mode in a special way. We do not assume 
that the phase ip(t) is a small quantity, but allow it to 
undergo a Brownian-like motion under the action of the 
stochastic noise [46]. Substituting the above into Eq. (2) 
and going into Fourier space, the fluctuations S(j)kpo{k, t) 
separate from S4>k=o(t) and ip(t) in the linearized limit. 
One gets a set of equations 


d S4>k=o 


dt 

dip 

dt 


-2r Scj)k=o 


1 dWi^k=o 

Uri Jt 




2pL^/^ Scj)k=o 


dWii^k=o 

dt 


(19) 

( 20 ) 


where Nc{t) = + (5</'fc=oP and we split 

the noise into the real and imaginary part with 
{dWRj^k=odW^j j,^Q) = D^^dt. Equation (19) corre¬ 
sponds to the Ornstein-Uhlenbeck process while Eq. (20) 
describes phase diffusion. By stochastic integration of the 
first equation [47] we calculate the number fluctuations 
in the condensate 

An,.o = 2No^^^ = y (21) 

As usual for an interacting condensate, the number 
fluctuations scale with the square root of Nq. The 
mean number of particles in the condensate and out-of- 
condensate modes is given by (recall that Nq = l^opT'^) 


fiVo + 

D(p(p 

4r 

[ r 

4+r" 1 1 


for k = 0, 
for fc ^ 0, 


( 22 ) 


which should be compared with Eq. (16). The singularity 
at A: = 0 has been removed in Eq. (22), but an additional 
term /4r indicates that on average the noise slightly 
increases the number of particles in the condensate. 

For consistency, the size of the box must be suffi¬ 
ciently small so that the fraction of particles in the 
k = 0 mode is considerable. The number of parti¬ 
cles in out-of-condensate modes can be estimated by re¬ 
placing the sum of over k by an integral over the 
modes from AA; to Ai^ax = 21 x1 L. From this estimate 


we get the self-consistency “small box” condition to be 
^ ^ ™ ID case and Llnit§^ < 

in the 2D case. 


V. NUMERICAL RESULTS 


In this subsection we investigate in detail the effects 
of the breakdown of the adiabatic approximation on the 
momentum distribution and spatial correlation functions 
of a one-dimensional condensate. To this end we inves¬ 
tigate numerically the steady states of the system using 
both the ID version of the SGPE model (2) and the open- 
dissipative Gross-Pitaevskii model (4)-(5) with the fluc¬ 
tuations included. We compare these numerical results 
with the analytical predictions given by (16), (17) and 
( 22 ). ^ ^ 

To investigate the breakdown of adiabaticity we con¬ 
sider various values of parameters for which the adia¬ 
baticity conditions (13) and (14) are either fulfilled or 
not. We clearly observe that the two models lead to the 
same results only in the adiabatic regime. We also find 
that the discrepancy between the SGPE and ODGPE 
models may be caused by large density fluctuations, 
which also invalidate the correspondence between the 
models. This behaviour is observed close to the limits 
of condensate stability (15) and P = Pth- 

Note that in the one-dimensional case the nonlin¬ 
ear coefficients scale with the confinement strength, 
eg. (P^°, 5 i°) = (P^°, gf^)/'J2ixd ?, if we assume a 
Gaussian transverse profile of |^/)|^ and ur of width d. 
In the case of a one-dimensional microwire [48], the pro¬ 
file width d is of the order of the microwire thickness. 
For convenience we also introduce dimensionless param¬ 
eters in the system (4)-(5) which allows us to describe 
the system using a limited set of parameters. By rescal¬ 
ing time, space, wave function amplitude and material 
coefficients as t = rt, x = xp = = 

i^P)~^nR = {^f3/T)R, (ff^°,g]f) = {KPlT){g,gR), 
{ic,1 r) = t- ^{7c,7r), P{x) = (1/^/3t)P(x), where 
^ = y/hT/2m*, while r and /3 are arbitrary scaling pa¬ 
rameters, we can rewrite the above equation in the di¬ 
mensionless form (from now on we omit the tildas for 
convenience) 


.dxp 
* dt 

duR 

dt 


d'^ i 

“WT + O (P'^R - ic) +3c|'0p + gR^R 

dx ^ 




P - {lR + R\'^P?)nR, 


(23) 


In particular, we may choose the time scaling r in such 
a way that 7 r = 1 without loss of generality. The norms 
of both fields = J and Nr = f URdx are 

scaled by the factor of /3, and the noise correlators be¬ 
come {dW {x)dW* {x')) = P^iPriR + 'yc)5x,x'- The 
parameters of the SGPE model (2) can be rescaled in an 
analogous way. 
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FIG. 1: (Color online). Stable and adiabatic regime. Left panel shows the position of the system on the stability-adiabaticity 
diagram. The black line corresponds to the stability condition (15) and the blue line to the adiabaticity condition (14). 
Middle panel shows momentum distribution in the steady state n“‘{k) = \ip{k)fAk for both models compared to the analytical 
prediction from (16), (22). On the right panel the calculated first-order correlation functions g^^\x) are displayed together with 
the analytical long-range fit (17). Parameters in dimensionless ODGPE units (23) are R — 1, gc = 0.4, gu = 2gc, P/Pth = 4, 
7c = 1, 7H = 1. /3 = 0.003, L = 600. 



FIG. 2: (Color online). Stable and non-adiabatic regime. Same as in Fig. 1, except for g — 20, gn = 2g, P/Pth = 2.5, 7 c = 0.5. 


The numerical simulations were performed on a grid 
of length L = 600 with spatial step size Aa; = 0.133. We 
solved the equation (2) and (4)-(5) with the correspond¬ 
ing parameters given by the conditions (10). We com¬ 
pared the calculated momentum distribution n®'*(fc) and 
the first-order correlation function (x) in the steady 
state after long time evolution, starting from the initial 
state with (l){x,t = 0) = (po and nB.{x,t = 0) = 'yclR- 
We found that typically the result does not depend on 
the initial conditions. However, for certain values of pa¬ 
rameters, close to the stability threshold (15), an empty 
initial state (j>{x, 0) = 0 led to a qualitatively different 
evolution (see the next subsection for details). 

A. Stable-adiabatic regime 

In Figure 1 we show the results in the stable-adiabatic 
case, i.e. for parameters fulfilling both the stability (15) 
and adiabaticity conditions (13)-(14). The left panel 
shows the position on the phase diagram. The middle 
panel displays the numerically obtained momentum dis¬ 
tribution n®®(fc) = |'!/'(fc)pAfc from 15 simulations of both 
models, compared with the analytical prediction of the 


momentum “tail” {k > 0) of (16) and (22). While the 
numerical box size L is so large that we are in a quasi¬ 
condensate rather than condensate regime [25], and the 
fraction of particles at fc = 0 is very small, the momentum 
tail distribution closely follows the analytical predictions. 
This follows from the fact that in the quasicondensate 
regime the Bogoliubov approximation can still be applied 
to a slice of the system where the quasicondensate phase 
is approximately constant. 

The discrepancy between the SGPE and ODGPE 
models is visible only at high momenta, where the 
first (momentum-dependent) adiabaticity condition (12) 
breaks down. Indeed, this condition gives the upper limit 
/c <C 8 , which coincides with the value of k at which the 
ODGPE results start to deviate from the SGPE and an¬ 
alytical predictions. We note that this limit corresponds 
to very high momentum values. 

On the other hand, the g^'^\x) function follows the 
analytically predicted long-range trend (17) very closely 
for both models. The difference between the two is in the 
different value of the constant in front of the exponent 
in (17), which can be again attributed to the breakdown 
of adiabaticity at very high momenta, related to short 
distances. 
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FIG. 3: (Color online). Steady state close to the condensation threshold, exhibiting large fluctuations. Same as in Fig. 1, except 
g = 2.5, OR = 2g, P/Pth = 1.05, yc = 0.1. 
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FIG. 4: (Color online). Density evolution corresponding to 
Fig. 3, revealing dark solitons in the dynamical steady state 


of the ODGPE. 


Further, we performed a series of numerical tests in 
the stable-adiabatic regime with other values of param¬ 
eters. With increasing ratio ^c/lR, moving towards the 
stability threshold (15), the differences between analyti¬ 
cal and numerical results become more notable. Beyond 
the stability threshold, which corresponds to go = 0, the 
comparison is no longer possible since the analytical for¬ 
mulas are valid only for the case of repulsive effective 
interactions, go > 0. 


B. Stable-nonadiabatic regime 

Figure 2 shows the results in the stable-nonadiabatic 
case, when the adiabaticity conditions (13) and (14) are 
not fulfilled. In this case we increased the strength of the 
interactions gc and qr. There is a visible discrepancy 
between the models both in the momentum distribution 
and the first-order correlation function, which marks the 
breakdown of the adiabatic approximation. Note that 
the results are presented in the logarithmic scale, and 


the actual difference between the calculated averages of 
n^^{k) differ signihcantly. The SGPE result still follows 
the Bogoliubov analytical prediction closely. The corre¬ 
lation function only slightly differs from the analytical 
prediction for both models, as shown in the right panel. 


C. Large fluctuations 


Importantly, as shown in previous sections, the adia¬ 
baticity conditions (12)-(14) are valid only under the as¬ 
sumption that the system is in a steady state with small 
fluctuations. Moreover, without this assumption the re¬ 
duction of the generalized CGLE (9) to the GGLE or 
SGPE models is not possible. We illustrate the situation 
in which the fluctuations are large in Fig. 3, where the 
pumping was chosen slightly above the threshold value 
P ~ Pth- In this case even the generalized Bogoliubov 
approximations are not relevant and the analytical pre¬ 
dictions are incorrect. Moreover, the results provided 
by the SGPE and ODGPE models are qualitatively dif¬ 
ferent. The SGPE model predicts momentum distribu¬ 
tion n®'’(fc) that is practically independent of momentum 
and negligible spatial correlations. On the other hand, 
the ODGPE predicts a nontrivial momentum distribu¬ 
tion and decaying g^^\x) indicating an appearance of a 
finite spatial correlation length, not related to the Bo¬ 
goliubov prediction (17). 

We associate this spatial length scale with the sponta¬ 
neous appearance of dark structures depicted in Fig. 4. 
The figure shows density distribution in one randomly 
chosen realization of the truncated Wigner stochastic 
evolution (4)-(5), which can be interpreted as a single 
realization of the experiment [14, 31]. The structures 
appear to be related to dark solitons of the dissipative 
model [44, 49, 50]. We checked that each dark object cor¬ 
responds to an approximate tt phase jump of the phase 
of ipixjt). A detailed investigation of these excitations 
will be presented elsewhere. 

















FIG. 5: (Color online). System close to the stability threshold. In the right panel momentum distributions obtained from 
different initial conditions are shown (see text). Parameters are as in Fig. 2, except 7 c = 1. 


D. Stability threshold and the role of initial 
conditions 

In the majority part of the stability diagram, the ini¬ 
tial conditions given as a seed to the evolution do not 
influence the steady state properties. However, the sit¬ 
uation is drastically different for the ODGPE model in 
the vicinity of the stability threshold given by (15). An 
example is given in Fig. 5. Here, the momentum distri¬ 
bution is plotted for simulations starting from two differ¬ 
ent initial conditions, either ip{x,t = 0) = ipQ + ^(x) or 
'ip{x,t = 0) = C(x), where ^ represents a small uncorre¬ 
lated noise with a Gaussian distribution. In the first case, 
the system converges to a steady state that is very well 
described by the analytical Bogoliubov momentum distri¬ 
bution. In the second case, the system does not reach this 
state even after very long evolution, instead dwelling in 
a chaotic evolution with large density fluctuations. The 
“normal” and “chaotic” states are therefore metastable. 
This behavior can be generally observed in the vicinity 
of the stability limit, both on the stable and unstable 
side of the phase diagram of Fig. 5. We checked that this 
bistability persits even if a relaxation term (frequency- 
dependent pumping) is included the evolution equation 
(4), analogous as in [51]. Since we do not find any similar 
dynamics in the SGPF model, we conclude that it is also 
related to the breakdown of the adiabatic approximation. 
The investigation of these chaotic states will be a topic 
of a future study. 

VI. CONCLUSIONS 

In conclusion, we investigated the relation between 
the models commonly used in the literature to de¬ 
scribe nonresonantly pumped exciton-polariton conden¬ 
sates. The complex Ginzburg-Landau equation, and the 
equivalent (in the limit of small fluctuations) stochastic 


Gross-Pitaevskii equation were compared with the open- 
dissipative Gross-Pitaevskii equation which includes a 
separate equation for the reservoir density. The adia¬ 
batic approximation allows to reduce the latter to one of 
the single-equation models, under the assumptions that 
the condensate is close to the steady state and fluctu¬ 
ations are small. Additionally, three independent ana¬ 
lytical conditions for adiabaticity must be met simulta¬ 
neously. While spin degree of freedom was not taken 
into account in this work, the generalization to the spin- 
dependent case is straightforward. 

We investigated the corresponding stochastic models 
by comparing the numerical steady states with the ana¬ 
lytical predictions of the Bogoliubov approximation. We 
demonstrated how the zero-momentum singularity of the 
momentum distribution spectrum can be avoided by an 
appropriate generalization of the Bogoliubov approxima¬ 
tion. This also allowed for determination of the number 
fluctuations and condensate phase diffusion equation. 

The comparison of the models with and without a sep¬ 
arate equation for the reservoir demonstrated that close 
agreement between the two can be obtained only under 
the adiabatic conditions. Moreover, we showed that close 
to the limit of condensation or the limit of modulational 
stability, large fluctuations lead to qualitatively differ¬ 
ent results depending on the model used. These results 
show that special care must be taken when choosing the 
right model for describing exciton-polariton condensates 
in certain conditions. 
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